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ABSTRACT 


Day by day, the integration of decentralized generation in medium voltage 
networks becomes more important during the last years and even in the near 
future. This increase causes, at the same time, several negative effects and 
rarely positive impacts on the stability of the network. Therefore, this work 
aims at analyzing the impact of ambient temperature on radial distribution 
network parameters’ e.g.: voltage drop and stability voltage level (index). 
Based on MATLAB program, different analyses of distributed generation 
(DG) insertion influence’s on voltage drop in the radial distribution feeder, as 
well as the influence of climatic conditions such as ambient temperature on 
network parameters. The Integration of Photovoltaic DGs in MV networks 
can play an important role in reducing the global warming effect (in voltage 
drop, and voltage stability index) especially in radial distribution feeder. 
Furthermore, it protects network’s parameters if its location and power are 
well selected. 
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1. INTRODUCTION 

In the electricity sector the first step in understanding the impacts of climate change is based on 
knowing the certain terms such as increasing the ambient air temperature [1]. Around the world, 
environmental pollution has experienced a real acceleration in recent years; this acceleration has several 
forms such as the increase in heat of the gallop. In a desert region the heat can vary in normal cases between 
5 up to 50° C, and on this basis several systems are installed especially those affected by the ambient 
temperature such as electrical energy distribution systems. The data bases used by the operators in the tables 
of facilities are drawn up on the basis of 25° C and 50° C. But this reality will change over time as a result of 
global warming, for example; the ADRAR region (desert) recorded a maximum temperature on July 9, 2018 
which reached 65°C [2]. These new records recorded can seriously affect the stability of the electrical power 
distribution system in general, and the parameters dependent on the particular ambient temperature such as 
voltage drop. 

Therefore, the electrical distribution systems have to be re-studied, and get an overall idea of the 
increase in ambient temperature impacts on the voltage stability of the electrical system, and the use of DGs 
as solution which also have impacts on stability. In this context, several studies focused on the integration 
impacts of the decentralized generation in the distribution networks, these impacts can take several aspects. 
Some researchers were carried out on DG technology and its impacts on the energy system. These impacts 
may concern network planning, losses, voltage, energy quality and reliability [3-8]. The large number of DG 
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penetration has an influence on voltage disruption, voltage drop in distribution systems particularly with self- 
disconnection of loads [9-10]. Therefore, when DGs are integrated, some standards should be respected such 
as the insertion place as well as the DGs number and its size. The location and technology type of DGs 
(active and reactive energy injection, active energy injection only, active energy injection and absorbs 
reactive energy) effects the energy losses [11-16]. Some researches determine the optimal penetration level to 
30% to minimize the DGs intermittent effects such as wind turbines and photovoltaics [17-18]. With 
the great growth of the installation of the DGs, others focus on algorithms and models that make it possible 
managing and controlling the distribution system and analyze the loads flow such as active power and 
reactive communication systems [19-20]. Several studies are done on voltage stability index and its relation 
to DG type, to rate contribution and to DGs integrated parameters [21-29]. 

Some studies show the influence of the ambient temperature on the electric system, as the impact 
study of the increase of ambient temperature and harmonics on the transformers which reduces its life time 
[30]. Moreover, the impact of the climatic changes (increase of ambient temperature) on energy system 
infrastructures can reduce the DGs capacity especially in the desert regions as Western United States [31-32]. 
Other studies show the changes in resistance values of the transmission lines caused by the temperature 
ambient are non-negligible effects [33]. This paper focuses on the ambient temperature impact on voltage 
stability of the radial distribution systems (voltage stability index, voltage drop) of electric energy system in 
the desert regions, as well as the advantage of the insertion of distributed generation. For This, this paper is 
divided into two main parts; the first consists of concepts and equations that relate the drop and stability of 
the voltage at ambient temperature. As for the second part it depends on a case study on a 4bus feeder test 
system and IEEE 33 bus radial distribution system. 


2. PROBLEM FORMULATION 
In this paper, the objective functions of ambient temperature to reconfiguration problem are voltage 
drop and voltage stability index which are calculated for each case in radial distributed systems. 


2.1. Voltage drop (VD) 

Voltage drops (VD) is the difference between the voltage at the source (or sending end) bus (Vs), 
which is assumed to be a fixed value, and the voltage at the load (or receiving end) bus (VL), which varies as 
a function of the load current (l1). [34] 


2.2. Factors affecting voltage drop 

One can classify the factors affecting the voltage drop in two categories. The first is dependent to 
load parameters that are; voltage, current and power factor, the second factor depends on line impedance 1.e.; 
length and resistance. This last is linked directly to the material nature (copper, aluminum ...), to 
the conductor cross-section and ambient temperature. 


2.3. Voltage drop calculations 

The determination of voltage drop value for lines inferior to 50 km in medium voltage network (U < 
40 kV) is based on two formulas: 

First formula that gives the exact result: 


VD = Vs + IRcos(@) + IXsin(@) — [Vs? — (IRcos(@) + [Xsin(6))?]*/? (1) 


Where: VD is the voltage drop, I is the line current (A), R is the circuit line resistance (Q), X: is the 
circuit line reactance (Q), and @ is the load power factor angle (degrees) Vs; source voltage, VL load voltage. 
The voltage drops VD that is calculated as a line-to-neutral voltage drop (one-way), and the line-to-line 
voltage drop. And a second approximate formula that nearly gives the same values 1s: 


VD = IRcos(@) + [Xsin(@) = IZ (2) 
Z=Rcos(@) + Xsin(@) (3) 


that is calculated for a single-phase system is 2VD; and the tine-to-tine voltage drop that is calculated for a 
three-phase system is V3VD. So, one can determine the percentage of the voltage drop by 
the following formula: [35] 


VD(%) = — * 100 (4) 
S 
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2.4. Ambient temperature impact on voltage drop 
The change in resistance as a function can be represented by the following formula: 


R2 = R1 [(k + T2)/(k + T1 )] (5) 


where: 
R1 is a conductor resistance in ohms (Q) at temperature T1 in degrees Celsius (°C) 
R2 is a conductor resistance in ohms (Q) at temperaturesT2 in degrees Celsius (°C) 
and k is the resistivity temperature coefficient at 0°C for the conductor; k= 234.5 for copper, Aluminum 228, 
lead 230, steal 202. [35] 
For T=T; et R&R; 
VD1= I*[Ri«cos (8) +X: *sin(0)] 
and for T=T2 , R= Ro 


R2 = Ry [(k + T,)/(k +T )] (6) 
So 

VD 2=I(R2«cos 8 +X +sin 8) (7) 
X; is independent of ambient temperature. Therefore, the (7) becomes; 

VD, = I[R,[(k + T,)/(k + T,)| * cos 0 + X, * sin] (8) 


and from (4) and (7) the percent voltage drop becomes 


k+T3 
k+T1 





VD.% = 1 [| R( )| «cos@ + X; *sin®] * 100/Vs (9) 


voltage drop equal V3VD for three phase system 
So, the percent voltage drop for three phase system formula becomes; 


k+T> 
k+T1 





VD,% = v3 1I [Ra ( )| «cos@ + X; *sin®] * 100/Vs (10) 


the (8) means that there is a relation between the ambient temperature and voltage drop variation. 


2.5. Voltage stability index (VSD 

Voltage stability is the ability of a power system to maintain steady acceptable voltages at all buses 
in the system under normal operating conditions and after being subjected to a disturbance, a system enters a 
state of voltage instability when a disturbance, increase in load demand, or change in system condition causes 
a progressive and uncontrollable drop in voltage [36]-[38]. 


2.6. Voltage stability index calculations 
Voltage stability index can determined in favor of this formula: 


VSI(r) = 2VéV?2 — Vi — 2V2(PR + QX) — |Z|? (P? + Q?) (11) 


when s is source, r is loads and R, X, (Z) impedance of line, P and Q is charges [39]. In reality there are other 
formula to calculate 


2.7. Ambient temperature impact’s on voltage stability index 

The formula shows that voltage stability level has a big relation with the line impedance, voltage 
buses, and load power. The impedance value increases with increase of its resistance as the ambient 
temperature increases. This affect the voltage stability index in the different buses of the power system (see 
formula 5) So, formula 10 becomes: 


VSI(r) = 2VéV?2 — Vi —2V2(P*R1 [(k + T2)/(K+T1 )] + QX) —|Z'|?(P?7 +903 (11) 
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3. METHODOLOGY AND SIMULATION 

We will devise this part of simulation in tow devises, the first talk about the impact of ambient 
temperate in voltage drop value for case to IEEE 4-Bus Test Feeder, and who can limit its effects as shown in 
Table 1. The second simulation devise of take IEEE 33-Bus Test Feeder radial distribution system like 
example to effects of insertion of DGs on voltage stability index value, and selection of the good conditions 
to support voltage stability. 


Table 1. Parameters of IEEE 4 bus test feeder elements 
Source Line 1 (L1) Cable 1 (L2) Transformer Load 
Ss=100 MVA Ra&= 0.09 OĒQ/km Ræ&= 0.09 OQ/km 12,47/4,16 Kv  Vı=4,16 kV 
Vs= 12,47 kV Xa= 0.52 O/km X= 0.52 Q/km = Rey= 0.132 Q l=1200 A 
Rs= 0.12 Q Imax= 2000 A Imax= 2000 A Xc= 0.620 X/R =15 
Xs= 0.48 Q L,=3 km L2= 0,672 km Pf = 90% 


3.1. Fist case: IEEE 4-bus test feeder 


To carry out a simulation of the distribution network during hot days, the 8th of June temperature 
model will be applied (Figure 1) on the model IEEE 4 Node test feeder (Figure 2). 
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Figure 1. Temperature in ADRAR region for the 9" of Juley 2018 
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Figure 2. IEEE 4-bus test feeder 


By using MATLAB program, the variations of voltage drop in the 4 buses during 24 hours can be 
determined (see Figure 3). 

This figure shows that the ambient temperature variance influences the impedance of the overhead 
lines and finally over the entire distribution network. But there are two cases; a positive case that is attached 
to the ambient temperature decrease, and a second considered negative that is the case of the ambient 
temperature increase, it is considered by the growth of the percentage (%) of the voltage drop. 
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Figure 3. rate of voltage drop change in bus1,2,3,4 during 24h 


3.2. Advantage of DG in voltage drop reduction 

In order to see the nature and the intensity of the temperature influence on voltage drop calculations 
and especially in the future medium voltage networks where connected DGs may be found. Based on that, 
the following application in two different seasons is taken. The first one in winter where the change in 
temperature from -5 °C up to 20°C, the second one in summer that varies from 20°C up to 60°C, and 
sometimes exceeds 65°C in some Sahara regions and sub-of the Sahara. Consequently, the study interval 
should be taken from -5°C up to 70°C. In the first case, an application is carried out on a simple network that 
contains a source, a transmission line MT and a load, the voltage drop for the bus 2 is therefore calculated 
(see Figure 4 (a)). In the second case, a DG source is connected to see how the changes made influence (see 
Figure 4 (b)). The DG contributes in the load supply with 50%. And Table below presents data of 
the different elements of the network studied. Figure 5 determines the ambient temperature impact on 
the voltage drop value in bus 2. The calculations made by technicians are often made for an ambient 
temperature of 25 ° C where VD = 3.1%. It can be noted that the temperature can exceed this threshold and 
the calculations made becomes very far from reality. The insertion of a DG source at the proximity of 
the bus 2 plays a very positive role in the backup of the VD in these standards (VD is ~ 7% for MV [40]). 


Contribution of DGs in the stability and voltage drop reduction for... (Bouafia Abdelkader) 


982 o ISSN: 2088-8694 


Bus MV ii Bus2 


Oa] 


Source 
(a) 


Bus MV T Bus2 


Onm 
Source 
(b) (26) 


Figure 4. Simple medium voltage network with and without DG connection 
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Figure 5. Voltage drop change in bus 2 under ambient temperature influence with and without DG insertion 


Table 2. Numerical data for simple circuit like figure 2 


Source DG Line 1 (L1) Load 
Ss=120 MVA Spo = 10 MVA Ry = 1,65 Q/km V,=11,5 kV 
Vs= 30 kV Vps= 30 kV X= 1,22 Q/km I,=105,7 A 
Rs= 0.12 Q Rpc= 0,80 Q Imax= 2000 A X/R =15 
Xs= 0.48 Q Xpc= 0,32 Q L,\=3 km Pf=0,95 


3.3. Second case: IEEE 33-bus radial distribution 

In order to confirm the impact of DG insertion and its location on the voltage stability index (level) 
under ambient temperature, an MV 33 radial bus distribution network (Figure 6) is used. This distribution 
system contains 33 bus (nodes) and a 3.715 MVA, 2.3 MVAR, 12.47kV power source. For further data see 
the Appendix (A & B). 

The measured voltage in the different bus bars under ambient temperatures is: -5°C, 25°C and 65°C 
respectively. These temperatures represent the two ends of T2 and 25 °C represents T1. The results obtained 
are presented in Figure 7. The figure shows the degradation of the voltage level along the line, especially in 
buses 18 and 33 under the effect of the ambient temperature increase. This also causes a decrease in the level 
of voltage stability index (see Figure 8). The future MV distribution networks are networks with high 
possibility of the DGs insertion in one or more points. This type of networks suffers from climate heating but 
has the advantage of being in areas of very high annual average solar radiation, so it is very advantageous to 
install DGs of photovoltaic PV types. In fact, it has firstly been dealt with this type of problem by trying 
integrating a PV type DG with a 0.5MW active power in different IEEE 33 bus bars radial distribution 
system, the results obtained are shown in Figure 9. The tests indicate that the two sites that are really 
improving are inserts 18 and 33. 
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Figure 6. Single line diagram of IEEE 33-bus radial distribution 
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Figure 7. Voltage drop in differences buses for -5, 25, and 70°C respectively 
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Figure 8. Voltage stability index in buses for different ambient temperatures 
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Figure 9. Result to insertion of DG-PV in different buses 
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3.4. DG advantage in voltage stability 

In order to understand the impact of DG insertion or in a MV distribution system, the previous 
example in the two main points which give good performance in reducing the voltage drop will be used. Four 
photovoltaic type DGs with different sizes were chosen, see Table 3. These DGs supposed to generate active 
power only. 


Table 3. Distributed generation type and parametrers 


G Active Power (MVA) Reactive Power (MVAR) 
DG1-PV 0.5 - 
DG2-PV 1.0 - 
DG3-PV 2.0 - 
DG4-PV 4.0 - 


4. RESULTS AND DISCUSSION 

Thirteen (13) possible scenarios of the DGs insertion in IEEE 33-bus radial distribution see table 4. 
At the beginning the minimum and the maximum of the voltage stability index (SVI) for two different 
ambient temperatures 25 °C and 70 °C is measured, then the divergence min-max (SVI) for each scenario 
will be calculated. The results obtained are shown in Figure 10. 

Three categories of scenarios can be distinguished; better, average and bad. The 10 and 11 scenarios 
are classified as best category that are related to the insertion of two small DGs (0.5 & 1MW) in the bus 18 
and 33, the max of this category does not exceed 1 (SVI <1 stability), and the divergence min-max is 
minimum as well as the impact of the ambient temperature appeared almost zero. for the average category 
one can classify the scenarios 2,3,4 and 5 which are related to the insertion of a single small DG either in 
the buses 18 or 33, the max of SVI and almost 1, and the divergence Min-Max is under the initial value 
(without DG), the impact of the ambient temperature is minimum. concerning the bad category it is found 
the scenarios related to the insertion of one to the two DGs of large sizes (2 or 4 MW), which are 
the scenarios (6, 7, 8, 9, 12 and 13), the max usually exceeds 1 (SVI > 1 instability) and the divergence Min- 
Max is greater than the initial value, and the impact of the ambient temperature is remarkable. From these 
results it can be noted that the insertion of the small DGs at the end of the branches plays a significant role in 
the stability of the distribution networks and limits the impact of the ambient temperature on the installation. 


Table 4. Presents min and max of Stability voltage index for differencts scenarios 


Scenario DG used Position Min SVI for 25°C Min SVI Max SVI Max SVI 
for 70°C for 25°C for 70°C 

Scenario 1 Without - 0,66116286 0,62277933 1 1 

Scenario 2 DGI1-PV bus 18 0,72400454 0,70386996 1 1 

Scenario 3 DGI1-PV bus 33 0,68629247 0,65581334 1 1 

Scenario 4 DG2-PV bus 18 0,74888752 0,72897958 1 1 

Scenario 5 DG2-PV bus 33 0,74888752 0,68456904 1 l 

Scenario 6 DG3-PV bus 18 0,79633411 0,78088314 1,26719789 1,30687082 

Scenario 7 DG3-PV bus 33 0,76190067 0,74142979 1,0306693 1,05094967 

Scenario 8 DG4-PV bus 18 0,87798095 0,87212255 2,00982316 2,16617314 

Scenario 9 DG4-PV bus 33 0,85830433 0,84886693  1,41000027 1,48224003 

Scenario 10 Tow DG1-PV  busl8 & 33  0,80079868 0,78886449 1l 1 

Scenario | 1 Tow DG2-PV _ bus18 & 33 0,89593788 0,89625162 1,00037338  1,01569951 

Scenario 12 Tow DG3-PV bus18 & 33 0,94841805 0,94400568  1,40370693 1,47522629 

Scenario 13 Tow DG4-PV _ bus18 & 33 0,96031167 0,95452289 2,37014846 2,6318921 
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Figure 10. divergence min-max of Stability voltage index to differencts scenarios for 25°C and 70°C 


5. CONCLUSION 

The electric energy distribution system is usually exposed to external impacts such as the ambient 
temperature that affects the internal stability parameters of the electric system that is usually voltage. For 
future electric distribution networks, the climate change effects should be taken into consideration during 
the phase of installation and especially for the calculations of the transmission lines impedances by keeping 
the ambient temperature superior or equal to 75°C. The insertion of the DGs remains a good means to fight 
against the impact of these temperature peaks that often touches the Sahara zones. The DGs insertion in a 
radial distribution network and its impact on voltage drop and on voltage drop indices under the effect of 
the difficult climatic conditions (extreme hot temperature zones) has been studied in this paper. 

It has been shown that DGs can play a very important role in the improvement of the voltage 
stability index, as well as limitation of the voltage drop for future electric distribution networks, selection of 
the size, number and the insertion points of the DGs that well reduces of voltage drop and give a good 
stability. The impact of the ambient temperature on voltage drop (VD) and on voltage stability indication 
(VSI) that appears minimal in the different chosen scenarios will become more important in the most 
complex electric distribution networks and especially during peak days. According to the study, the optimal 
position where the DGs should be inserted in the radial distribution networks is at the ends of the branches as 
the 33 buses case been studied. The latter deals with the insertion of two photovoltaic type DGs with medium 
or small dimensions. This has shown a good efficiency in terms of stability for radial distribution network 
and minimized considerably the increase of the ambient temperature impact. 
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Appendix: 
Table A. Load data of the IEEE-33 radial distribution Table B. Line data of the IEEE 33-bus radial 
system distribution system 
Load Location RealLoad Reactive Load (kVAR) Branch From To Length Impedance 
(bus) (kW) bus bus (km) Resistance _Reactance 
Load2 2 100 60 l 1 2 1 0.0922 0.0470 
Load3 3 90 40 2 2 3 1 0.4930 0.2511 
Load4 4 120 80 3 3 4 1 0.3660 0.1864 
Load5 5 60 30 4 4 5 1 0.3811 0.1941 
Load6 6 60 20 5 5 6 1 0.8190 0.7070 
Load7 7 200 100 6 6 7 1 0.1872 0.6188 
Load& 8 200 100 7 7 8 1 1.7114 1.2351 
Load9 9 60 20 8 8 9 1 1.0300 0.7400 
Load10 10 60 20 9 9 10 1 1.0440 0.7400 
Load1 1 11 45 30 10 10 11 1 0.1966 0.0650 
Load12 12 60 35 11 11 12 1 0.3744 0.1238 
Load13 13 60 35 12 12 13 1 1.4680 1.1550 
Load14 14 120 80 13 13 14 1 0.5416 0.7129 
Load15 15 60 10 14 14 15 1 0.5910 0.5260 
Load16 16 60 20 15 15 16 1 0.7463 0.5450 
Load17 17 60 20 16 16 17 1 1.2890 1.7210 
Load18 18 90 40 17 17 18 1 0.7320 0.5740 
Load19 19 90 40 18 2 19 1 0.1640 0.1565 
Load20 20 90 40 19 19 20 1 1.5042 1.3554 
Load21 21 90 40 20 20 21 1 0.4095 0.4784 
Load22 22 90 40 21 21 22 1 0.7089 0.9373 
Load23 23 90 50 22 3 23 1 0.4512 0.3083 
Load24 24 420 200 23 23 24 1 0.8980 0.7091 
Load25 25 420 200 24 24 25 1 0.8960 0.7011 
Load26 26 60 25 25 6 26 1 0.2030 0.1034 
Load27 27 60 25 26 26 27 1 0.2842 0.1447 
Load28 28 60 20 27 21 28 1 1.0590 0.9337 
Load29 29 120 70 28 28 29 1 0.8042 0.7006 
Load30 30 200 600 29 29 30 1 0.5075 0.2585 
Load31 31 150 70 30 30 31 1 0.9744 0.9630 
Load32 32 210 100 31 31 32 1 0.3105 0.3619 
Load33 33 60 40 32 32 33 l 0.3410 0.5302 
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